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as  a  tool  in  reservoir  water  quality  management. 
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INITIAL  EVALUATION  OF  A  COMPUTER  ALGORITHM  FOR  SIMULATING 
ANAEROBIC  CONDITIONS  IN  RESERVOIRS  AND  LAKES 


PART  I :  INTRODUCTION 


Study  Objective 


1.  The  potential  for  development  of  low  dissolved  oxygen  (DO)  or 

anoxic  conditions  is  a  major  concern  for  reservoir  managers.  Low  DO 

concentrations  in  impoundments  can  lead  to  the  generation  of  undesirable 

chemical  products  of  anaerobic  biogeochemical  processes.  These  prod- 

+2  +2 

ucts,  which  include  dissolved  manganese  (Mn  ),  iron  (Fe  ),  and  sulfur 
-2 

(S  ),  along  with  the  low  DO  condition,  can  degrade  recreational  aesthe¬ 
tics  and  fisheries  and  can  cause  taste,  odor,  and  staining  problems  that 
increase  water  treatment  costs.  Hydrogen  sulfide  (H^S)  releases  are 
known  to  cause  illness  In  reservoir  project  operators,  and  also  accele¬ 
rate  corrosion  of  outlet  works  (Keeley  et  al.  1978).  In  addition,  in- 

-3  + 

creased  regeneration  of  phosphorus  (PO^  -P)  and  ammonia  (NH^-N)  may 
stimulate  nuisance  algal  blooms. 

2.  In  order  to  simulate  effects  of  low  oxygen  conditions  on  water 
quality  in  planned  and  operational  reservoirs,  the  Corps  of  Engineers 
(CE)  has  incorporated  new  algorithms,  in  the  form  of  13  subroutines, 
into  its  one-dimensional  numerical  model  of  water  quality,  CE-QUAL-R1 
(Environmental  Laboratory  1982).  The  objective  of  this  report  is  to 
document  the  development  and  testing  of  the  new  subroutines  in  simulat¬ 
ing  anaerobic  processes  as  they  pertain  to  reservoirs.  This  work  was 
performed  during  early  stages  of  model  development  and  does  not  reflect 
recent  improvements,  some  of  which  came  in  response  to  results  of  this 
study . 


Model  Background 
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3.  CE-QUAL-R1  describes  the  vertical  distribution  of  thermal 


energy  and  biological  and  chemical  material  for  preimpoundment  and  post- 
impoundment  studies  of  water  quality.  Conceptually,  CE-QUAL-Rl  depicts 
a  reservoir  as  a  series  of  horizontal  layers  (Figure  1)  within  each  of 
which  thermal  energy  and  mass  are  uniformly  distributed.  Mathemati¬ 
cally,  the  model's  structure  is  based  on  variable  horizontal  layer 
thickness,  which  depends  on  Inflowing  and  outflowing  water  volume; 
inflowing  water  is  distributed  among  the  layers  by  density.  Vertical 
transport  of  energy  and  mass  is  achieved  through  entrainment  and  turbu¬ 
lent  diffusion. 

4.  The  model  can  simulate  the  dynamics  of  many  water  quality  var¬ 
iables  including  temperature,  DO,  phytoplankton,  zooplankton,  dissolved 
solids,  pH,  alkalinity,  dissolved  inorganic  carbon,  suspended  solids, 
dissolved  organic  matter,  phosphorus,  nitrogen,  iron,  manganese,  and 
sulfur.  Water  quality  simulations  are  performed  for  each  layer  of  the 
model.  A  complete  description  is  given  in  the  CE-QUAL-Rl  User's  Manual 
(Environmental  Laboratory  1982). 

Innovations  for  Simulating  Anaerobic  Conditions 


5.  An  overriding  consideration  in  developing  the  new,  anaerobic 
algorithms  is  to  keep  the  model  as  simple  as  possible  while  adequately 
representing  the  system.  To  this  end,  the  chemistry  has  been  simpli¬ 
fied,  although  the  number  of  additional  variables  may  seem  to  indicate 
otherwise. 

6.  The  original  concepts  for  subroutines  to  deal  with  anaerobic 
processes  within  the  framework  of  CE-QUAL-Rl  appear  in  previous  tech¬ 
nical  reports  (Brannon  et  al.  1978,  Gunnison  and  Brannon  1981).  Those 
reports  also  contain  valuable  background  information  and  data.  Addi¬ 
tional  important  concepts  for  the  new  algorithms  are  based  on  data  from 
Environmental  and  Water  Quality  Operational  Studies  (FWQOS)  field 
studies  and  on  the  model's  one-dimensional  assumption. 

7.  Two  simplifying  conceptual  innovations  mark  the  operation  of 
the  anaerobic  algorithms.  The  first  is  the  use  of  a  DO  concentration  in 
the  water  column  as  a  cue  to  initiate  anaerobic  processes.  Intensive 
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scrutiny  of  field  data  from  several  locations  suggested  that  a  reason¬ 
able  correlation  existed  between  oxygen  decline  in  the  water  column  and 
the  release  and  reduction  of  anaerobic  materials.  Lum,  Armstrong,  and 
Fruh  (1981)  and  Fillos  and  Molof  (1972)  have  reached  similar  conclu¬ 


sions.  Based  on  field  data  from  DeGray  Lake  in  the  1979  model  calibra¬ 
tion  year,  threshold  DO  concentration  selected  was  0.5  mg*i  .  That  is, 
examination  of  field  data  indicates  that  when  declining  DO  concentra¬ 
tions  fall  to  0.5  mg*i  dissolved  anaerobic  materials  appear  in  the 
water  samples.  This  threshold  value  subsumes  many  organic  decay  or 
oxygen-depleting  processes  which  may  vary  among  different  systems. 

8.  It  may  seem  that  redox  chemistry  axioms  are  being  violated 
deliberately  by  invoking  the  DO  cue.  However,  the  primary  redox  pro¬ 
cesses  of  concern  in  these  subroutines  are  those  resulting  in  sediment 
release  of  dissolved  materials.  Since  we  have  no  data  on  redox  state  or 
dynamic  chemical  conditions  in  the  sediments,  DO  in  the  water  column  is 
used  as  a  correlate  of  sediment  redox  condition.  The  assumption  that 
this  cue  is  valid  in  the  water  column,  as  well,  may  be  questioned. 

9.  Subsequent  to  the  threshold  DO  concentration  being  reached, 

releases  of  so-called  "anaerobic  materials"  (PO.^-P,  NH^-N ,  Mn+^,  Fe+^, 
-2  4  4 

and  S  )  from  sediments  and  reduction  processes  in  the  water  column  do 

not  occur  sequentially  as  they  might  if  redox  potentials  were  used  as 
cues.  Thus,  the  second  major  innovation  consists  of  starting  the  anaer¬ 
obic  processes  of  sediment  release  and  chemical  reduction  simultane¬ 
ously.  Each  process,  however,  has  its  own  rate  which  is  fixed  initially 
but  may  be  affected  by  temperature. 

10.  The  model  predicts  water  quality  as  though  a  specific  point 
in  a  given  layer  represents  the  whole  layer.  A  real  system  would  likely 
be  more  diverse  in  composition;  that  is,  spatially  separated  points  at 
the  same  elevation  are  unlikely  to  be  chemically  identical.  Allowing 
simultaneous  release  from  the  sediments  of  the  various  anaerobic  mate¬ 
rials  at  different  rates  is  intended  to  account  for  the  true  heterogene¬ 
ity  of  the  reservoir  system  within  the  context  of  a  one-dimensional 
model.  Initiating  anaerobic  processes  simultaneously  and  calibrating 
rates  to  simulate  field  conditions  may  resolve  electrochemical  parodoxes 


such  as  the  contemporaneous  appearance  of  dissolved  oxygen  and  dissolved 
reduced  metals. 


Anaerobic  Subroutines 


11.  Thirteen  new  mathematical  subroutines  in  CE-QUAL-R1  calculate 
water  quality  variables  specifically  pertinent  to  anoxic  or  low  DO  con¬ 
ditions.  Each  subroutine  deals  with  an  individual  chemical  component  of 
water  quality  and  its  physical,  chemical,  and  biological  interactions 
with  other  water  quality  constituents  in  the  same  or  adjacent  layers. 

12.  In  a  conceptualized  horizontal  layer,  anaerobic  constituents 
(or  compartments)  occur  either  in  the  water  column  or  in  the  sediments 
(Figure  2).  In  the  water  column,  anaerobic  water  quality  constituents 
are  generally  present  either  in  the  oxidized  (particulate)  or  reduced 

(dissolved)  state,  depending  on  DO  concentration  and  duration  of  anoxia; 

_2 

oxidized,  dissolved  sulfur  (SO^  )  and  reduced,  particulate  iron  sulfide 
are  exceptions. 

13.  Anaerobic  materials  in  the  sediments  comprise  6  of  the  13 
compartments.  These  sediment  compartments  have  two  primary  functions. 
First,  they  serve  as  a  source  for  most  of  the  anaerobic  materials. 
Second,  they  provide  a  method  of  accounting  for  net  fluxes  of  anaerobic 
materials  between  the  sediments  and  overlying  waters.  Within  the  sedi¬ 
ments,  actual  oxidation  states  are  immaterial  for  simulation  purposes 
since  sediment  electrochemistry  is  not  explicitly  considered.  Again, 
the  important  idea  is  that  anaerobic  materials  are  released  simultane¬ 
ously  when  DO  concentrations  fall  below  the  threshold  value. 

14.  The  following  brief  summary  describes  the  events  depicted  in 
Figure  2  as  DO  falls  below  the  threshold  concentration. 

a.  Ammonia  (NH+-N) ,  dissolved  (reduced)  manganese  (Mn  ^) , 

~  4  +2-2 

dissolved  (reduced)  iron  (Fe  ),  sulfide  (S  ),  and 

_3 

dissolved  phosphate  (PO^  -P)  are  released  from  the 
sediments  at  specific  rates. 
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b.  Particulate  manganese  (Mn(IV)),  particulate  iron 

_2 

(Fe(III)),  and  sulfate  (SO^  )  may  be  reduced  in  the  water 
column. 

c.  Dissolved  iron  and  sulfide  may  form  insoluble  iron  sul¬ 
fide,  which  may  precipitate  through  the  water  column  to 
the  sediments. 

15.  To  ensure  that  the  supply  of  anaerobic  materials  from  the 
sediment  is  not  depleted  prematurely,  a  coefficient  defining  an  active 
sediment  depth  is  included.  When  this  value  is  multiplied  by  concentra¬ 
tion  in  the  sediment  and  the  sediment  surface  area,  the  product  is  the 
available  mass  of  a  given  sediment  component  in  a  particular  layer. 

This  allowed  a  check  on  mass  balance. 

16.  When  DO  is  replenished  and  exceeds  the  threshold  value,  oxi¬ 
dation  processes  predominate:  ammonia  oxidizes  to  nitrite  and  then  to 
nitrate;  particulate  manganese  and  iron  are  generated  from  dissolved; 
sulfide  is  oxidized  to  sulfate;  and  iron  sulfide  is  oxidized  to  sulfate 
and  particulate  iron. 

17.  The  subroutines  representing  anaerobic  processes  contain  both 
stoichiometric  and  kinetic  coefficients.  The  stoichiometric  coeffi¬ 
cients  determine  oxygen  depletion  resulting  from  oxidation  processes. 
These  reactions  and  their  stoichiometric  oxygen  requirements  (grams 
consumed :gram  material  oxidized)  are  listed  in  Table  1  by  the  name  as¬ 
signed  each  coefficient  in  the  computer  program. 

18.  In  calculating  the  dynamics  of  a  given  compartment  describing 
an  anaerobic  material,  additional  rate  coefficients  are  required  for 
oxidation-reduction  processes,  sediment  release  of  reduced  substances, 
and  settling  velocities  of  particulates.  Reaction  rate  coefficients 
express  the  fraction  of  the  oxidation  or  reduction  reaction  completed 
per  day.  Release  rates  refer  to  the  amount  of  anaerobic  material 
released  across  a  square  meter  of  sediment  surface  per  day.  A  settling 
velocity,  in  meters/day,  describes  the  downward  movement  of  a  precipi¬ 
tated  chemical  compound  (oxidized  iron  or  manganese,  or  iron  sulfide)  in 
the  water  column.  Other  simplifications  involve  the  use  of  zero-order 


Table  1 

Stoichiometric  Coefficients  Affecting  Dissolved 
Oxygen  Concentration 


kinetics  to  describe  rates  of  release  of  anaerobic  materials  from  sedi¬ 
ments  and  first-order  reaction  kinetics  for  oxidation-reduction  trans¬ 
formations. 

19.  Reaction  rates  are  supplied  initially  by  the  CE-QUAL-R1  user, 
but  may  be  affected  by  temperature.  Because  most  chemical  and  biologi¬ 
cal  reaction  rates  are  temperature  dependent,  a  mathematical  function  is 
incorporated  in  the  model  to  determine  the  effect  of  this  relationship 
(Thornton  and  Lessem  1978).  Figure  3  shows  the  variations  of  the  rate 
multiplier  factors  from  user-supplied  temperatures  (T^  to  T^,)  and  multi¬ 
plier  coefficients  (K^  to  K^) ,  as  discussed  in  the  RMULT  section  of  the 
CF.-0UAL-R1  User's  Manual  (Environmental  Laboratory  1982). 

20.  As  an  example  of  the  calculations  involved,  the  following 

equation  describes  the  various  processes  affecting  the  concentration  of 
+2 

Fe  in  a  given  layer  in  the  model's  formulation: 

Tt  (VC>  •  k  (DA  i  M)  +  EQinCl„  -  «outC  +  Wr.dVC* 

-  WoxTC  -  Y3WC  +  V6Kr.lA. 

where 

3 

V  «  layer  volume  in  m 

+2  -3 

C  -  Fe  concentration  in  g*m 

+2  2  -1 
D  «*  Fe  diffusion  coefficient  in  m  •hr 

A  ■  layer  area 

AZ  =  layer  thickness  in  m 

3  -1 

Q.  ■  layer  inflow  in  m  *hr 

ln  3  -1 

Q  =  layer  outflow  in  m  *hr 
out 

“  temperature  rate  multipliers  (dimensionless) 

K  «  reduction  rate  of  Fe  (III),  hr  * 
red 

C*  *  Fe  (III)  concentration  in  g*m 

+2  -1 

K  **  oxidation  rate  of  Fe  ,  hr 

ox  +2  -2 

K_  -  Fe  reaction  rate  with  S  to  form  FeS 

u  +2-2-1 
K  ,  »  sediment  release  rate  of  Fe  in  g*m  *hr 

rel  +2  2 
Ag  ■*  area  of  sediment  releasing  Fe  ,  m 


Thus,  the  first  right-hand  term  in  the  equation  describes  changes  in 
+2 

Fe  due  to  diffusion  between  layers;  the  second  term  accounts  for 

additions  from  inflowing  waters  summed  over  all  tributaries;  the  next 

+2 

term  subtracts  losses  through  layer  outflow;  the  fourth  term  adds  Fe 

+2 

derived  from  Fe(III)  reduction;  the  next  term  determines  losses  of  Fe 

via  oxidation  to  Fe(III);  the  sixth  term  calculates  the  amount  of  dis- 
+2  -2 

solved  Fe  lost  through  reaction  with  S  to  form  FeS;  and  the  final 
+2 

term  increments  Fe  from  sediment  release.  The  CE-QUAL-R1  User's 
Manual  contains  complete  descriptions  of  the  chemical  equations  and  the 
details  of  their  programming. 


PART  II:  CALIBRATION  OF  ANAEROBIC  SUBROUTINES 


Site  Characteristics 


21.  In  1972,  DeGray  Lake  was  selected  as  a  field  study  site  for 
research  on  water  quality  modeling  under  the  CE  Environmental  Impact 
Research  Program  (EIRP) .  The  lake  was  chosen  because  of  the  availabil¬ 
ity  of  preimpoundment  and  postimpoundment  baseline  data,  its  multiple 
outlet  structure,  and  pumped  storage  hydropower  capability.  In  1977, 
the  EWQOS  program  selected  DeGray  Lake  as  one  of  four  reservoir  study 
sites  to  provide,  among  other  objectives,  information  needed  in  develop¬ 
ing  and  verifying  water  quality  predictive  techniques.  In  this  study, 
data  from  1979  were  used  to  calibrate  the  model,  and  data  from  1980  were 
used  to  confirm  the  model. 

Description 

22.  DeGray  Lake,  created  in  1970  by  the  Impoundment  of  the  Caddo 
River,  is  located  near  Arkadelphia  in  south-central  Arkansas.  Its 
watershed  is  approximately  30  percent  agricultural,  with  the  rest  for¬ 
ested.  DeGray  Lake  is  a  multipurpose  reservoir  with  flood  control, 
power  generation,  recreation,  water  supply,  and  low-flow  augmentation  of 
the  Ouachita  River  as  project  purposes.  The  lake  is  a  highly  dendritic 
system,  approximately  32  km  long,  with  a  mean  depth  of  15  m  and  a  max¬ 
imum  depth  of  60  m.  The  project  has  pumped  storage  capability  and  can 
withdraw  water  at  three  levels  in  the  lake  (Ford  and  Stein  1982) .  Aver¬ 
age  hydraulic  residence  time  is  0.8  year  (Thornton  et  al.  1982). 
Climatology 

23.  DeGray  Lake  experiences  a  modified  continental  climate.  Sum¬ 
mers  tend  to  be  long,  warm,  and  humid  with  air  temperatures  occasionally 
exceeding  40°  C.  Short,  mild  winters  are  occasionally  affected  by  cold 
polar  fronts  capable  of  dropping  temperatures  to  -20°  C.  Mean  annual 
temperature  is  around  17°  C. 

24.  Precipitation  is  relatively  evenly  distributed  throughout  the 
year,  with  the  March-May  period  accounting  for  about  one-third  of  the 
total.  Summer  rains  generally  consist  of  scattered  showers  and 


thunderstorms  that  may  be  locally  intense.  In  summer  and  fall,  there 
may  be  extended  periods  of  drought  or  near-drought  conditions.  Late 
fall  is  usually  marked  by  increasing  rainfall.  Average  annual  precipi¬ 
tation  at  Arkadelphia  is  134  cm. 

25.  Late  winter  and  early  spring  are  marked  by  the  highest  aver¬ 
age  daily  wind  speed,  while  lowest  wind  speeds  are  experienced  in  mid- 
to  late-summer. 

26.  The  meteorological  background  information  given  above  was 
based  on  the  report  of  Ford  and  Stein  (1984).  The  following  material 
comes  from  examination  of  dry  bulb  temperatures,  fractional  cloud  cover, 
and  wind  speed  data  taken  at  Little  Rock,  Ark. ,  approximately  80  km 
northeast.  From  these  data,  some  general  comparisons  can  be  made  be¬ 
tween  the  calibration  (1979)  and  verification  (1980)  years. 

27.  Winter  1979  was  considerably  colder  than  1980.  In  late  win¬ 
ter  1979,  warming  occurred  with  an  abrupt  jump  of  about  5-6°  C,  before 
continuing  a  steady,  gradual  rise,  barely  reaching  30°  C  in  summer.  In 
1980,  the  gradual  warming  trend  was  not  marked  by  an  initial  fast  rise 
but  reached  peaks  in  excess  of  30°  C  several  times  during  the  summer. 
Autumn  temperatures  declined  somewhat  faster  in  1980  than  1979. 

28.  Cloud  cover  may  reflect  the  temperature  regimes  for  1979  and 
1980.  Cloud  cover  seems  to  have  been  slightly  less  in  early  1980  than 
1979.  Summer  1980  had  significantly  less  cloud  cover  resulting  in 
greater  insolation  and  higher  temperatures.  The  latter  part  of  1980  was 
perhaps  slightly  cloudier  than  1979. 

29.  Average  daily  wind  speed,  a  significant  factor  in  reservoir 
mixing,  averaged  about  12-14  km*hr  *  during  the  first  half  of  1979, 
dropped  to  about  8-10  km* hr  during  the  summer,  and  picked  up  again  to 
about  12-13  km*hr  *  for  the  rest  of  the  year;  highest  winds  rarely  ex¬ 
ceeded  20  km*hr  In  1980,  wind  speeds  were  about  15  km*hr  *  during 
the  first  third  of  the  year,  dropped  to  about  12  km*hr  *  in  the  middle 
of  the  year,  and  rose  to  around  13-14  km* hr  *  during  the  final  third. 
Highest  wind  speeds  in  1980  reached  or  surpassed  25  km*hr  at  the  be¬ 
ginning  and  end  of  the  year.  Overall,  wind  speeds  were  probably  a  bit 
higher  in  1980  than  1979. 


30.  These  trends  or  differences  are  relatively  crude  comparisons, 
but  they  do  serve  to  point  out  that  year-to-year  weather  variations  may 
be  significant.  The  ultimate  effect  of  these  differences  on  ecological 
simulations  is  at  times  difficult  to  assess. 

Stratif ication-destratif ication  cycle 

31.  DeGray  Lake  is  a  warm,  monomictic  system  that  typically 
stratifies  by  late  March.  Following  thermal  stratification,  metalimne- 
tic  oxygen  minima  develop  in  the  lake's  upper  reaches  and  eventually 
extend  to  the  rest  of  the  reservoir.  Oxygen  depletion  occurs  rapidly 
and  may  result  in  total  hypolimnetic  anoxia  in  the  relatively  shallow, 
up-reservoir  regions.  The  main  body  of  the  reservoir  does  not  exhibit 
such  intense  anoxia.  This  general  stratification  behavior  typifies  many 
CE  reservoirs. 

32.  Historically,  complete  vertical  mixing  first  occurs  in  the 
upper  part  of  the  reservoir  in  October  and  in  the  main  body  by  December; 
however,  overturn  expected  in  late  fall  1979  was  delayed  until  February 
1980.  Hydrometeorological  conditions  left  the  system  poised  for  mixing 
(isothermal  but  chemically  stratified)  but  lacking  an  energy  impulse 
(e.g.,  wind)  to  initiate  mixing. 


33.  Fortnightly  sampling  at  several  locations  in  the  DeGray  Lake 
system  provided  data  necessary  for  computer  simulations.  Other  inten¬ 
sive,  short-term  studies  of  storm  events,  primary  productivity,  and 
areal  variation  of  selected  water  quality  constituents  complemented  the 
routine  sampling  program.  Water  samples  were  taken  and  observations 
made  at  locations  representing  the  major  inflow  and  outflow  sites  in  the 
system  as  well  as  at  several  in-pool  locations. 

Water  quality  data 

34.  Water  quality  sampling  in  1979  was  performed  on  a  2-week 
basis,  with  the  exception  of  the  last  quarter  of  the  year  when  weekly 


samples  were  taken  (for  sampling  techniques,  see  Kennedy,  Montgomery, 
and  James  1983) .  In  order  to  make  the  best  possible  use  of  available 
data,  weekly  updates  were  used  for  simulation  Input,  and  missing 
alternate  weekly  values  for  the  first  9  months  were  determined  by  linear 
Interpolation.  Missing  analytical  data  points  were  also  determined  in 
this  manner.  In  order  to  account  for  several  minor  sources  of  inflowing 
water  affecting  water  balance,  a  second  tributary  was  included,  with 
water  quality  characteristics  identical  to  the  major  inflow,  the  Caddo 
River . 

35.  Several  of  the  model's  variables  reflect  theoretical  con¬ 
siderations.  Their  values  and  those  of  their  associated  coefficients 
are  either  theoretical  or  based  on  assumptions  concerning,  and  calcula¬ 
tions  based  on,  other  data.  Variables  that  are  calculated  from  field 
data  Include  detritus,  dissolved  organic  matter  (DOM),  and  suspended 
solids.  Anaerobic  constituents  of  sediments — nitrogen,  iron,  manganese, 
phosphate,  and  sulfur — are  defined  as  variables  and  given  values  calcu¬ 
lated  from  bulk  sediment  analysis;  these  variables  are  included  mainly 
to  track  fluxes  between  water  and  sediments. 

Calibration  Overview 

36.  It  is  of  the  utmost  importance  to  predict  accurately  the 
temporal  and  spatial  occurrence  of  the  threshold  DO  concentration  that 
cues  the  initiation  of  anaerobic  processes.  Photosynthesis,  algal  res¬ 
piration,  and  biologically  mediated  decay  processes,  as  well  as  physical 
processes  that  affect  mixing,  may  significantly  influence  simulated 
oxygen  profiles.  In  particular,  parameters  associated  with  algal  physi¬ 
ology  and  growth  and  organic  decay  appear  critical  in  calibrating  bio¬ 
logical  effects  on  oxygen  concentrations  in  the  water  column,  as 
indicated  in  preliminary  work. 

37.  The  influence  on  the  anaerobic  subroutines  of  the  biological 
processes  associated  with  phytoplankton  are  beyond  the  scope  of  this 
report . 


Oxygen  Calibration 

38.  Once  the  model  has  been  calibrated  to  provide  satisfactory 
thermal  profiles,  the  graphics  utility  of  CE-QUAL-R1  should  be  used  to 
examine  DO  profiles  for  the  time  period  under  consideration.  Experience 
and  intuition  may  play  important  roles  in  deciding  which  coefficients 
and  variables  to  examine  first  to  achieve  the  best  data  fit.  Table  2 
lists  initial  coefficient  values  used  in  this  calibration  study.  Values 
for  coefficients  not  originating  with  the  anaerobic  algorithms  were  used 
in  earlier  work  and  may  be  found  in  the  User's  Manual  (Environmental 

Table  2 

Initial  Values  of  Coefficients  Used  in  Calibration  Study  of 
Anaerobic  Processes  in  CE-QUAL-R1 

_ Coefficient 

Oxygen  cue 

Sediment  thickness 

Mn(IV)  settling  rate 

Mn(IV)  reduction  rate 
+2 

Mn  sediment  release  rate 
+2 

Mn  oxidation  rate 


Initial  Value 
0.5  mg •£  * 

5.0  cm 
0.05  m*day  * 

0.02  day-1 

-2  -1 
0. 10  g*m  *day 

0.00  day-1 


-1 


Laboratory  1982).  FeS  reactions  were  not  examined  in  this  report,  so 
associated  initial  values  are  set  to  0.0. 

39.  Reservoirs  frequently  develop  metalimnetic  DO  minima.  As 
thermal  stratification  develops  and  intensifies,  anoxic  conditions  may 
first  occur  in  the  metalimnion  and  then  extend  more  or  less  gradually  to 
the  hypo limn ion,  which  has  been  "sealed  off"  by  the  metalimnetic  density 
gradient.  Therefore,  when  these  conditions  are  known  to  occur,  the  ini¬ 
tial  step  in  oxygen  calibration  should  concern  accurately  predicting  the 
onset  of  anoxia  in  the  metalimnion.  Table  3  summarizes  the  simulations 
performed  in  this  study. 

40.  Figure  4  demonstrates  the  development  of  the  DO  minimum  in 
DeGray  Lake  simulations  and  field  observations.  The  simulation  compares 
favorably  with  the  observed  data  with  regard  to  the  initiation  of  anoxic 
conditions  in  early  August.  The  model  slightly  underestimates  the 
amount  of  oxygen  present  in  the  hypolimnion  after  stratification  is  well 
established  in  the  summer.  Changes  in  the  decay  rate  of  an  organic  com¬ 
ponent,  such  as  DOM,  detritus,  or  sediment,  would  seem  to  provide  a  rem¬ 
edy  for  this  discrepancy.  However,  to  decrease  overall  rates  would  also 
affect  epilimnetic  oxygen  concentrations,  which  appear  satisfactory. 
Therefore,  modification  of  the  temperature  rate  multiplier  (see  Fig¬ 
ure  3)  for  DOM  to  decrease  decay  rates  in  the  cooler  hypolimnetlc  waters 
should  consequently  result  in  higher  oxygen  concentrations. 
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Table  3 

Summary  of  Simulation  Figures 


Type  of  Simulation  or 

Graph  of : 

Coefficient  Tested 

Value  Used 

°2 

Initial  simulation 

Values  in  Table  2 

»2 

DOM  rate  multiplier 

0.10 

°2 

DOM  rate  multiplier 

0.04 

0.0040  day  * 

°2 

DOM  decay  rate 

°2 

Sediment  decay  rate 

0.009  day"1 

0.10  day"1 

°2 

Detritus  decay  rate 

NH.-N 

4 

Initial  simulation 

-2  -1 

0.01  g  N»m  *day 

nhI-n 

4 

NH*!"-N  release  rate 

4 

-2  -1 

0.00  g  N*m  *day 

PO43-1 

P 

Initial  simulation 

-2  -1 

0.001  g  P*m  *day 

-3 

PO.  -1 
4 

P 

-3 

PO^  -P  release  rate 

-2  -1 

0.000  g  P*m  »day 

Mn+2 

Initial  simulation 

-2  -1 

0.1  g  Mn*m  *day 

Total 

Mn 

Initial  simulation 

-2  -1 

0.1  g  Mn*m  *day 

Mn+2 

Mn+2  release  rate 

-2  -1 

0.05  g  Mn*m  »day 

Mn+2 

Mn+2  release  rate 

-2  -1 

0.20  g  Mn*m  »day 

Total 

Mn 

Mn+2  release  rate 

-2  -1 

0.05  g  Mn*m  »day 

Total 

Mn 

w  +2 

Mn  release  rate 

-2  -1 

0.20  g  Mn*m  »day 

Mn+2 

Mn(IV)  reduction  rate 

0.00  day  1 

Total 

Mn 

Mn(IV)  reduction  rate 

0.00  day  1 

Mn+2 

Mn(IV)  reduction  rate 

0.02  day  1 

Total 

Mn 

Mn(IV)  reduction  rate 

0.02  day  1 

„  +2 
Fe 

Initial  simulation 

-2  -1 

0.10  g  Fe *m  »day 

Total 

Fe 

Initial  simulation 

-2  -1 

0. 10  g  Fe*m  «day 

v  +2 

le 

+2 

Fe  release  rate 

-2  -1 

0.05  g  Fe»m  *day 

Total 

Fe 

_ ,  +2  , 

Fe  release  rate 

-2  -1 

0.05  g  Fe«m  *day 

_  +2 
Fe 

t,  +2 

Fe  release  rate 

-2  -1 

0.20  g  Fe*m  *day 

Total 

Fe 

„  +2 

Fe  release  rate 

-2  -1 

0.20  g  Fe*m  *day 

(Continued) 


Table  3  (Concluded) 


Figure 

No. 


Graph  of: 


Type  of  Simulation  or 
Coefficient  Tested 


Value  Used 


Fe^ 

Fe(III)  reduction  rate 

0.00  day 

Total 

Fe 

Fe(III)  reduction  rate 

0.00  day  * 

Fe+2 

Fe(III)  reduction  rate 

0.06  day  * 

Total 

Fe 

Fe(III)  reduction  rate 

0.06  day  * 

s'2 

Initial  simulation 

-5  -2 

10  g  Smb  •< 

s'2 

S  release  rate 

-4  -2 

10  g  S«m  •' 

s“2 

-2 

SO.  reduction  rate 

0.000  day-1 

s-2 

-2 

SO.  reduction  rate 

0.002  day-1 

s-2 

-2 

SO^  reduction  rate 

0.000  day-1 

-2 

S  release  rate 

0.1  g  S*m  ‘’•d. 

°2 

1980  simulation 

as  in  Table  2 

Mn+2 

1980  simulation 

as  in  Table  2 

Fe+2 

1980  simulation 

as  in  Table  2 

s'2 

1980  simulation 

as  in  Table  2 

°2 

Anaerobic  versus 

as  in  Table  2 

nonanaerobic 

4  SEPTEMBER 


2  OCTOBER 


30  OCTOBER 


41.  For  example,  the  lower  DOM  temperature  rate  multiplier  is 
initially  set  equal  to  0.12.  As  a  first  test,  lowering  the  value  to 
0.10  should  provide  an  indication  of  the  system's  sensitivity  to  this 
parameter.  Figure  5  presents  selected  simulated  and  observed  oxygen 
profiles.  Differences  from  the  original  profiles  (Figure  4)  are  mar¬ 
ginal.  More  drastic  changes,  e.g.  setting  the  temperature  rate  multi¬ 
plier  to  0.04,  yield  variable  results  (Figure  6).  The  simulation  curve 
fits  the  observed  profile  very  well  just  below  the  oxygen  minimum;  how¬ 
ever,  in  the  hypolimnion,  the  fit  is  considerably  worse.  In  this  case, 
then,  it  seems  best  to  leave  the  rate  at  0.12  and  to  attempt  to  improve 
the  simulation  by  altering  other  parameters. 

42.  Another  constant  that  may  have  an  impact  on  the  rate  of 
oxygen  depletion  is  the  DOM  decay  rate,  initially  set  at  0.0045  day  ^ . 
Reducing  this  rate  to  0.0040  day  ^  significantly  slows  the  dissolved 
oxygen  loss  (Figure  7).  While  this  change  does  generate  a  good  fit 
below  the  metalimnetic  oxygen  minimum,  it  also  delays  the  onset  of 
anoxia  in  the  metalimnion  and  causes  the  hypolimnion  to  experience 
higher  than  observed  oxygen  levels  at  the  end  of  the  year.  Results  such 
as  these  are  somewhat  ambiguous;  since  the  0.0045  day  1  value  gives  the 
better  approach  to  anoxia  initiation,  it  will  be  used  even  though  a 
value  of  0.0040  day  ^  may  generally  improve  the  hypolimnetic  oxygen  sim¬ 
ulation.  In  any  event,  it  is  obvious  that  DO  is  very  sensitive  to 
manipulations  of  this  rate. 

43.  Organic  sediment,  too,  comprises  a  potentially  large  oxygen 
sink  through  its  decay.  Sediment  is  one  of  the  model's  variables  for 
which  field  data  are  not  used  because  of  the  lack  of  information  on  the 
depth  of  sediment  interaction  with  the  water  column.  The  organic  sedi¬ 
ment  value  is  set  sufficiently  large  to  prevent  depletion.  The  sediment 
decay  rate  thus  may  exert  a  sizable  influence  on  the  rate  of  oxygen 
depletion.  Figure  8  shows  the  effect  of  a  50-percent  increase  in  the 
rate,  from  0.006  to  0.009  day  ^ .  Through  the  onset  of  metalimnetic 
anoxia,  the  effect  appears  minimal,  primarily  causing  a  slight  increase 
in  oxygen  loss  in  the  hypolimnion  which  persists  through  the  calendar 
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Figure  5.  Selected  DeGray  Lake  oxygen  profiles  in  1979  with  the 
DOM  temperature  rate  multiplier  set.  at  0,  10  (P  =»  predicted  value 

0  *  observed) 
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Figure  7.  Selected  DeGray  Lake  oxygen  profiles  in  1979  with  the 
DOM  temperature  rate  multiplier  set  at  0.004  (P  ■  predicted  value 

0  ■  observed) 
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44.  Similarly,  increases  in  the  rate  of  detritus  decay  from  0.006 
to  0.10  day  have  an  insignificant  effect  on  oxygen  profiles 

(Figure  9). 

45.  In  short,  for  this  particular  data  set,  manipulations  of  DOM 
decay  most  strongly  affect  DO  concentrations  when  compared  with  other 
organic  material  decay  rates.  While  conditions  may  vary  among  data 
sets,  this  suggests,  at  minimum,  that  DOM  decay  should  be  the  first  rate 
to  test  when  calibrating  other  data  sets.  Finer  adjustments  may  be  pos¬ 
sible  with  detritus  and  sediment  decay  rates. 
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Figure  9.  Selected  DeGray  Lake  oxygen  profiles  in  1979  with  the 

detritus  decay  rate  set  at  0.10  day  (P  **  predicted  value;  0  » 

observed) 


Calibration  of  Anaerobic  Materials 

Ammonium  and  orthophosphate 

46.  The  anaerobic  subroutines  deal  with  the  effects  of  anoxic 

conditions  on  five  elements.  Two  of  these  elements,  nitrogen  (in  the 

+  -3 

form  of  ammonia,  NH^-N)  and  phosphorus  (as  orthophosphate,  PO^  -P) ,  are 

also  Included  as  variables  in  the  "aerobic”  portions  of  CE-QUAL-R1; 
decay  and  other  aerobic  processes  may  act  to  mask  the  effects  of  anaero¬ 
bic  release. 

47.  Comparison  between  output  (Figure  10)  from  an  initial 


simulation  with  an  anaerobic  NH^-N  release  rate  of  0.01  g  N*m  *day 
and  a  second  run  with  a  zero  release  rate  (Figure  11)  reveals  two 
important  results.  First,  the  model  overestimates  NH+-N  concentrations 

throughout  the  year  irrespective  of  anaerobic  conditions.  Second,  the 

+  -2-1 
decrease  in  NH^-N  concentration  in  the  0.00  g  N*m  »day  release  rate 

simulation  demonstrates  that  a  slight  change  in  release  rate  may 

significantly  affect  concentrations  in  those  layers  experiencing 

anaerobic  conditions;  however,  the  potential  impact  of  this  latter 

observation  is  uncertain  when  one  considers  the  overall  high  simulated 

NH^-N  concentrations. 

4  -3 

48.  P0^  -P  simulations  (Figure  12)  also  overestimate  total  con¬ 
centrations  (except  at  the  surface,  where  comparison  is  quite  good). 

Even  when  the  phosphorus  release  rate  from  sediments  under  anaerobic 
conditions  is  set  to  zero,  improvement  is  minimal  (Figure  13). 

49.  The  discrepancies  between  predicted  and  observed  values  for 

+  -3 

NH^-N  and  PO^  -P  cloud  this  aspect  of  the  evaluation  of  the  anaerobic 
subroutines;  because  these  vital  nutrients  are  available  in  excess,  sub¬ 
sequent,  interrelated  phenomena  of  algal  growth,  photosynthesis,  and 
respiration  may  be  distorted.  It  may  be  that  the  instantaneous  disper¬ 
sal  into  the  water  column  of  inflowing  nutrients,  implicit  in  the  one¬ 
dimensional  assumption,  is  at  fault.  Nutrient  adsorption  to  suspended 
solids  may  also  contribute  to  these  discrepancies.  Program  modifica¬ 
tions  to  appropriately  model  the  nutrients’  dispersal  to  better  reflect 
actual  reservoir  conditions  by  adding  nutrient  adsorption  and  settling 
algorithms  to  the  model  were  incorporated  following  completion  of  this 
study  (Wlosinski  and  Collins  1985) . 

Manganese 

50.  In  the  simulated  reservoir  system,  manganese  may  occur  as  the 

+2 

dissolved,  reduced  metal  (referred  to  as  Mn  )  or  as  total  manganese 


(the  sum  of  dissolved  and  particulate  manganese,  assumed  to  be  the  oxy- 

hydroxide,  Mn(IV)OOH).  Following  the  onset  of  anoxia  and  during  strati- 

+2 

fication,  sediment  release  and  Mn  reduction  rates  control  manganese 
concentrations . 

51.  After  metalimnetic  anoxic  conditions  develop,  computer 
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Figure  11.  Selected  NH.-N  profiles  for  DeGray  Lake  in  1979  with 
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sediment  release  rate  =  0.00  g  N*m  ‘day  (P  =  predicted  value; 

0  *  observed) 
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Figure  14.  Selected  dissolved  Mn  profiles  for  PeGray  Lake 
in  1979  (P  ■  predicted  value;  0  **  observed) 
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Figure  15.  Selected  total  Mn  profiles  for  DeGray  Lake 
in  1979  (P  ■  predicted  value;  0  ■  observed) 


predictions  match  measured  values  quite  well  in  the  anoxic  zone  (Fig- 

+2 

ures  14  and  15).  Maximum  simulated  Mn  concentrations  occur  in  late 
November  on  a  date  when  measured  values  are  incongruously  low  (Fig¬ 
ure  14g);  after  that  dip,  concentrations  return  to  higher  levels 

(Figure  14h) .  Total  Mn  concentrations  develop  similarly,  peaking  in 

+2 

mid-December  when  some  of  the  dissolved  Mn  begins  to  oxidize  (Fig¬ 
ure  15e).  Observed  hypolimnetic  values  of  total  and  dissolved  manganese 
are  higher  than  in  simulations.  Observed  DO  concentrations  are  about 

1-2  mg*Jl  Standard  redox  chemistry  principles  would  deny  the  presence 
+2 

of  Mn  at  the  observed  hypolimnetic  oxygen  concentrations,  which  more 

+2  -1 

than  suffice  to  prevent  sediment  release  of  Mn  with  the  0.5  mg  D0*i. 
threshold. 

+2 

52.  The  nominal  value  for  Mn  release  from  sediments  is  0.1  g 
-2-1 

Mn*m  *day  .  For  heuristic  purposes,  several  figures  are  presented  in 
which  this  value  is  halved  (Figures  16  and  17)  or  doubled  (Figures  18 
and  19) .  Concentrations  of  dissolved  and  total  manganese  correspond¬ 
ingly  decrease  by  one-half  or  double.  No  effect  is  seen  in  the 
hypo limn ion. 

53.  Figures  20-23  depict  effects  of  varying  Mn(IV)  reduction 

rates.  Setting  the  rate  to  0.00  day  has  no  apparent  influence  on  dis- 
+2 

solved  Mn  (Figure  20)  and  only  a  minor  effect  on  total  manganese  in 
the  metalimnion  (Figure  21).  Doubling  the  initial  reduction  rate  of 
0.02  day  *  changes  neither  dissolved  nor  total  manganese  concentrations 
(Figures  22  and  23,  respectively). 

54.  In  short,  the  manganese  sediment  release  rate  appears  to  be 
the  principal  controlling  factor  affecting  manganese  concentration  in 
the  model. 

Iron 

55.  Iron  chemistry,  with  the  potential  for  reaction  between 
reduced  iron  and  sulfide,  presents  a  mere  complex  picture  than  manga¬ 
nese.  As  in  the  case  of  manganese,  the  model  adequately  simulates  dis¬ 
solved  and  total  iron  in  the  anoxic  metalimnion  (Figures  24  and  25) . 

Like  manganese,  however,  the  problem  of  high  observed  hypolimnetic  iron 
concentrations  recurs  in  spite  of  DO  levels  in  excess  of  0.5  mg  Oj*! 
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Figure  17.  Selected  dissolved  Mn  profiles  with  sediment  release 
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rate  «  0.20  g  Mn  *m  *d  (P  -  predicted  value;  0  «  observed) 
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Figure  18.  Selected  total  Mn  profiles  with  sediment 
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release  rate  =  0.05  g  Mn  *m  *day  (P  «  predicted 
value;  0  =  observed) 
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Figure  19.  Selected  total  Mn  profiles 
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Selected  total  Mn  profiles  with  Mn  (IV)  reduction 
0.04  day  *  (P  *  predicted  value;  0  *  observed) 
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Figure  23.  Selected  total  Mn  profiles  with  Mn  (IV)  reduction 
rate  -  0.04  day  ^  (P  *  predicted  value;  0  *  observed) 
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Figure  25.  Selected  total  Fe  profiles  in  DeGray  Lake 
in  1979  (P  ■  predicted  value;  0  =  observed) 
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Figure  27.  Selected  total  Fe  profiles  with  sediment  release  rate 
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**  0.05  g  Fe  *m  *day  (P  ■  predicted  value;  0  ■  observed) 
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Figure  28.  Selected  dissolved  Fe  profiles  with  sediment  release 
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rate  ■  0.20  g  Fe  *m  *day  (P  ■  predicted  value;  0  ■  observed) 
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Figure  29.  Selected  total  Fe  profiles 
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release  rate  ■  0.20  g  Fe  *in  *day  (P  1 
value;  0  m  observed) 
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There  may  exist  a  thin  anoxic  zone  along  the  sediment-water  Interface 

that  could  not  be  detected  during  sampling.  Halving  or  doubling  the 
+2 

Fe  release  rate  affects  dissolved  and  total  iron  in  direct  proportion 
to  the  rate  change  (Figures  26-29) . 

56.  Lowering  the  iron  reduction  rate  to  0.00  day  *  apparently  has 
no  effect  on  dissolved  or  total  iron  concentrations  (Figures  30  and  31). 
Doubling  the  rate  to  0.06  day  *  also  does  not  alter  the  iron  profiles 
(Figures  32  and  33) . 

57.  Again,  like  manganese,  iron  concentrations  appear  most  sig¬ 
nificantly  affected  by  release  rates  from  the  sediments. 

Sulfur 

_2 

58.  Sulfur,  in  the  reduced  form  of  sulfide  (S  ) ,  is  primarily 

controlled  by  its  release  rate  from  sediment,  reaction  with  reduced  iron 

_2 

to  form  insoluble  iron  sulfide,  and  SO,  reduction.  Figure  34  depicts 

4  _5  _2  -2  -1 

sulfide  profiles  for  a  release  rate  of  10  g  S  *m  *day  .  Consider¬ 
ing  the  wide  variability  in  observed  data,  the  simulations  seem  quite 
reasonable.  Simulated  sulfide  levels  peak  in  late  November  (Figure  34d) 

and  drop  to  nearly  zero  by  late  December  (Figure  34f ) ,  when  the  upper 

-2 

water  layers  have  become  oxygenated;  observed  S  levels  drop  in  Novem¬ 
ber  (Figure  34d)  and  Increase  again  in  December.  As  in  the  iron  and 
manganese  profiles,  the  model  does  not  account  for  the  presence  of  some 
hypolimnetic  sulfide. 

59.  Increasing  the  sulfide  sediment  release  rate  from  10  g 

-2  -2  -1  -4 

S  *m  *day  to  10  g  has  no  effect  on  sulfide  profiles,  implying  that 

the  model  is  quite  insensitive  to  this  parameter  (Figure  35) . 

-2  -1 

60.  Setting  the  S0^  reduction  rate  to  0.000  day  effectively 

eliminates  all  sulfide  from  the  water  column  (Figure  36).  Doubling  the 
initial  reduction  rate  from  0.001  to  0.002  day  doubles  the  sulfide 
concentration  (Figure  37). 

_2 

61.  Thus,  it  would  seem  that  S0^  reduction  dominates  sulfide 

formation  dynamics,  with  sediment  release  an  insignificant  factor  in  the 

-2  -1 
model.  However,  if  the  SO,  reduction  rate  is  set  to  0.000  day  and 

4  -2  -2  -1 

the  sulfide  release  rate  is  set  to  0.1  g  S  *m  *day  ,  sulfide  profiles 
of  significant  magnitude  can  be  generated  (Figure  38) .  Sulfate 
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reduction  initially  predominated,  largely  because  of  the  relatively  high 
concentration  of  sulfate  in  the  water  column. 

Discussion 

62.  This  discussion  attempts  to  examine  inadequacies  of  the  an¬ 
aerobic  subroutines,  to  assign  sources  of  error,  to  point  out  the  over¬ 
all  reliability  of  the  simulations,  and  to  suggest  where  future 
improvements  may  be  made.  For  the  most  part,  discrepancies  between 
predicted  and  observed  values  may  be  as  much  due  to  variability  in  sam¬ 
pling  and  analysis  and  limitations  of  the  one-dimensional  approach  as  to 
subroutine  formulation. 

+  -3 

63.  Release  of  NH^-N  and  PO^  -P  from  sediments,  generally  re¬ 
ferred  to  as  nutrient  regeneration,  should  constitute  an  extremely  im¬ 
portant  aspect  of  a  water  quality  model  simulating  anaerobic  processes. 

+  -3 

However,  too  high  NH^-N  and  PO^  -P  concentrations  under  simulated  aero¬ 
bic  conditions  preclude  observation  of  any  anaerobic  effects.  (Follow¬ 
ing  the  completion  of  this  report,  research  aimed  at  understanding  the 
mechanisms  transporting  and  recycling  nutrients  in  the  DeGray  Lake  sys¬ 
tem  has  resolved  these  discrepancies.) 

64.  Another  difficulty  arises  in  attempting  to  simulate  hypolim- 
netic  release  of  anaerobic  materials  with  the  0.5  mg  D0*£  *  cue.  Simu¬ 
lated  oxygen  concentrations  in  the  hypoliranion  fail  to  fall  low  enough 
to  initiate  sediment  release  in  the  deeper  hypoliranetic  layers.  Field 
data,  too,  would  indicate  that  DO  appears  too  high  to  allow  the  presence 
of  dissolved,  reduced  materials.  However,  all  anaerobic  materials  show 
up  in  spite  of  accepted  chemical  constraints.  Three  possible  explana¬ 
tions  for  the  occurrence  of  dissolved  metals  in  aerobic  hypolimnia  come 
to  mind.  First,  the  definition  of  dissolved  (=  reduced)  metal  is  an 
operational  one  based  on  passage  through  a  filter  of  specific  pore  size, 
not  on  electrochemical  considerations;  an  oxidized  molecule  may  be  able 
to  pass  through  the  filter.  Second,  the  metal  may  exist  in  a  dissolved 
state  in  association  with  an  organic  complex  or  colloid.  Third,  advec- 
tion  may  create  localized  pockets  of  reduced  materials  in  zones  of 


|  relatively  high  oxygen  concentration.  This  problem  is,  indeed,  a  dif- 

<  ficult  one,  which  limnologists  have  not  yet  resolved  and  which  may 

require  modification  of  the  model  to  allow  release  and  slow  oxidation  of 
dissolved  anaerobic  materials  in  aerobic  hypolimnia. 

65.  Two  aspects  of  the  field  data  contribute  to  less  than  ideal 
calibration  results.  First,  the  data  are  highly  variable  temporally  and 
spatially;  they  demonstrate  an  occasional  lack  of  continuity  in  trends. 
An  example  of  this  behavior  is  the  drop  in  concentration  of  anaerobic 
materials  in  late  November  (27  November  sample).  Since  the  decline  does 
not  appear  associated  with  any  increase  in  DO,  one  may  conjecture  that 
the  anomaly  resulted  from  sample  handling.  In  particular,  the  sulfide 
analysis  is  extremely  sensitive  to  the  slightest  oxygen  contamination. 

In  this  particular  case,  it  has  been  suggested  that  a  hydrological 
event,  missed  between  samples,  affected  the  trend.* 

66.  Second,  anaerobic  materials  occur  in  low  concentrations,  a 
factor  that  may  also  contribute  to  their  relatively  high  variability. 
Concentrations  of  Mn,  Fe,  and  S  are  reported  in  increments  of 
0.1  mg* i.  ,  although  analytical  techniques  should  yield  higher  degrees 
of  resolution.  Since  many  data  values  fall  near  zero  and  are  generally 
less  than  0.5  mg*Jl  1  for  any  anaerobic  material,  variability  of 
0.1-0. 2  mg*fc  *  in  field  observations  can  add  a  degree  of  uncertainty  to 
the  model  evaluation. 

67.  While  the  simulations  do  not  precisely  mirror  field  situa¬ 
tions,  they  do  match  observed  data  in  terms  of  initiation  of  release, 
magnitude  of  concentrations,  and  general  trends.  Agreement  between 
simulated  and  observed  conditions  is  actually  better  than  originally 
expected  when  taking  into  account  potential  sources  of  error  and  un¬ 
certainty,  such  as  simplicity  of  the  model  formulation,  the  one¬ 
dimensional  assumption,  use  of  estimated  parameters  and  coefficients, 
variability  in  environmental  sampling  conditions,  and  chemical  analyses. 


Personal  Communication,  1984,  J.  Nix,  Professor,  Ouachita  Baptist 
University,  Arkadelphia,  Ark. 


PART  III:  1980  CONFIRMATION  SIMULATION 


68.  The  same  coefficients  and  rate  constants  used  for  the  1979 
calibration  study  form  the  basis  for  a  simulation  of  1980  reservoir  con¬ 
ditions.  This  procedure,  referred  to  as  verification  or  confirmation, 
may  serve  several  purposes  and  its  results  may  be  Interpreted  in  differ¬ 
ent  ways.  In  general,  an  attempt  is  made  to  ascertain  the  validity  of 
the  model  examined,  or  of  any  of  its  parts,  including  coefficients.  In 
the  present  case,  the  1980  simulation  is  evaluated  to  determine  how  well 
the  new  anaerobic  subroutines  model  conditions  observed  in  DeGray  Lake 
for  a  year  other  than  that  used  in  the  calibration  study. 

69.  Since  DO  Integrates  physical,  chemical,  and  biological 
effects,  this  variable's  profiles  throw  light  on  the  adequacy  of  CE- 
QUAL-R1  to  model  the  reservoir  system.  Figure  39  traces  the  simulated 
development  of  anoxic  conditions  in  the  reservoir  and  compares  the  pre¬ 
dictions  with  observed  conditions.  As  the  graphs  indicate,  CE-QUAL-R1 
predicts  that  metalimnetic  anaerobic  conditions  (DO  s  0.5  mg*£  *)  begin 
in  late  July  (Figure  39g) ,  whereas  reservoir  DO  does  not  really  reach 
that  level  until  late  September  (Figure  39k). 

70.  Other  discrepancies  in  DO  profiles  exist.  Following  stable 
thermal  stratification  in  April,  epilimnetic  DO  predictions  match  ob¬ 
served  conditions  quite  well  until  July,  when  they  begin  to  diverge  sig¬ 
nificantly,  a  phenomenon  which  persists  and  becomes  more  marked  as  the 
year  progresses.  In  addition,  simulated  hypolimnetic  DO  never  drops  as 
low  as  field  data  indicate. 

71.  Physical  and  meteorological  conditions  may  combine  to  affect 
the  accuracy  of  DO  profiles  in  1980  simulations.  Wind  speed,  inflowing 
water  volume  from  tributaries,  and  thermal  structure  of  the  water  column 
differ  in  1979  and  1980;  these  three  environmental  factors  contribute  to 
the  calculation  of  the  diffusion  coefficient.  After  the  tapering  off  of 
high-volume  inflows  following  the  first  3-4  months  of  the  year  (1979  and 
1980),  wind  speed  is  the  primary  influence  on  the  diffusion  coefficient; 
the  average  wind  speed  in  1980  is  greater  than  in  1979,  possibly  result¬ 
ing  in  higher  diffusion  coefficients  in  1980.  Calibrating  the  diffusion 
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coefficient  functions  based  on  1979  data  may  lead  to  excessive  mixing 
and  increased  DO  throughout  the  water  column  in  1980.  (If  the  shelter¬ 
ing  coefficient  is  too  high,  similar  situations  may  develop.)  Another 
effect  of  increased  mixing  is  recharging  surface  waters  with  nutrients 
previously  depleted  by  algal  production.  Algal  productivity  Induced  by 
these  nutrient  additions  can  contribute  additional  0^  to  the  surface 
waters . 

+  -3 

72.  Since  NH^-N  and  PO^  -P  dynamics  could  not  be  calibrated 
properly  with  the  anaerobic  subroutines,  these  materials  will  not  be 
included  in  evaluating  1980  model  output. 

73.  Considering  that  simulated  anaerobic  conditions  start 

2  months  earlier  than  they  should,  and  assuming  all  other  factors  equal, 

one  would  expect  the  model  to  generate  anaerobic  materials  in  excess  of 

actual  concentrations.  This  expectation  is  borne  out.  Concentrations 
+2  +2 

of  dissolved  Mn  and  Fe  remain  higher  than  observed  values  until  some 

time  after  the  anoxia  criterion  (DO  ^  0.5  mg •£  *)  is  achieved  in  the 

reservoir  (Figures  40  and  41).  In  addition,  one  should  note  that  hypo- 

+2  +2 

limnetic  concentrations  of  Mn  and  Fe  remain  high  in  spite  of  high 

+2  +2 

DO.  This  situation  results  from  Mn  and  Fe  oxidation  rates  having 
been  set  to  0.00  day  *  in  the  1979  calibration  study.  In  the  present 
simulation,  initial  concentrations  of  these  anaerobic  materials  are  rel¬ 
atively  high,  but  oxidation  does  not  affect  them.  In  future  calibration 
work  on  different  systems,  it  may  be  advisable  to  maintain  at  least  a 

low,  but  positive,  default  value  for  oxidation  rates.  The  materials 
+2  +2 

Mn  and  Fe  experience  more  dynamic  changes  in  the  epllimnion  because 
mixing,  dilution,  and  outflow  affect  their  concentrations. 

74.  Sulfide,  on  the  other  hand,  does  have  a  nonzero  oxidation 

rate  and  so  does  not  appear  in  hypolimnetic  simulations.  As  in  the  1979 

_2 

calibration  comparisons,  observed  S  concentrations  vary  considerably 

with  no  discernible  trends  (Figure  42).  To  compound  the  difficulty  in 

studying  1980  simulations,  field  data  for  the  last  2  months  of  the  year 
-2 

do  not  include  S 

75.  In  spite  of  the  early  onset  of  anaerobic  conditions  predicted 

-1  +2  +2 
by  the  model  and  the  0.00  day  oxidation  rates  for  Mn  and  Fe  ,  the 
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Figure  41.  Selected  dissolved  Fe  profiles  in  1980 
(P  ■  predicted  value;  0  ■  observed) 
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Figure  42.  Selected  S  profiles  in  1980  (P  «*  predicted 
value;  0  »  observed) 


magnitudes  of  their  concentrations  never  become  unreasonably  high;  S 
concentrations,  likewise,  are  moderate.  Initially,  one  might  assume  the 
model  insensitive  to  the  increased  duration  and  extent  of  anaerobic  con¬ 
ditions.  More  likely,  the  calibrated  rates  are  extremely  low,  which  is 
in  line  with  observations  that  DeGray  Lake  at  the  sampling  site  never 
experiences  serious  anoxic  water  quality  problems.  Another  reason  for 
seemingly  low  accumulations  of  anaerobic  materials  in  the  metalimnion  is 
that  withdrawal  of  water  from  the  reservoir  is  largely  limited  to  the 
metalimnion  in  1980,  whereas  the  epilimnion  and  metalimnion  contributed 
in  1979.  Thus,  large  concentrations  of  anaerobic  materials  do  not  de¬ 
velop  in  the  metalimnion.  Futhermore,  water  budget  data  for  1974  to 
1980  (Ford  and  Stein  1984)  indicate  that  1979  was  an  atypical  year. 
Springtime  releases  were  considerably  greater  than  normal  (1980  was  a 
normal  year) ,  and  annual  release  volume  approximately  doubled  that  for 
1980.  In  spite  of  the  great  complexity  of  the  model  linking  physical, 
chemical,  and  biological  attributes  of  reservoirs,  the  simulations  yield 
results  which,  in  general,  reflect  actual  reservoir  conditions  and 
operations. 


69 


PART  IV:  IMPACT  OF  ANAEROBIC  SUBROUTINES  ON  CE-QUAL-R1 


76.  One  way  to  assess  a  new  version  of  a  model  is  to  compare  its 
predictions  with  its  predecessor's.  Oxygen,  the  only  variable  that  in¬ 
tegrates  effects  of  all  anaerobic  processes  in  the  model,  is  the  obvious 
choice  as  a  basis  for  comparison.  Since  the  earlier  version  of  the 
model  has  received  many  modifications  other  than  the  addition  of  anaero¬ 
bic  subroutines,  a  direct  comparison  of  oxygen  profiles  between  the  two 
versions  would  be  invalid;  any  number  of  changes  could  have  altered 
oxygen  profiles.  Therefore,  to  attempt  this  comparison,  only  the  most 
recent  version  of  CE-QUAL-R1  is  used,  either  with  or  without  implement¬ 
ing  the  anaerobic  subroutines.  Setting  the  value  of  the  DO  cue  to  -1.0 
effectively  bypasses  the  anaerobic  subroutines,  allowing  comparisons  be¬ 
tween  two  otherwise  identical  algorithms. 

77.  Figure  43  shows  how  similar  the  two  sets  of  oxygen  profiles 
are;  for  the  most  part,  they  are  identical.  However,  starting  on  2  Oc¬ 
tober  (Figure  43e)  and  continuing  through  the  year,  very  minor  differ¬ 
ences  appear  in  the  metallmnion  and  hypolimnion.  As  expected.  In  the 
simulations  which  incorporate  anaerobic  processes,  oxygen  concentrations 
are  slightly  lower  than  in  those  which  circumvent  anaerobic  processes, 
since  reduced  materials  create  an  extra  oxygen  demand.  Differences  are 
small  because  concentrations  of  reduced  materials  are  quite  low  and  may 
be  masked  by  oxygen  demands  of  organic  decay  processes.  Nevertheless, 
the  results  indicate  the  model  is  sensitive  enough  to  be  affected  by  the 
anaerobic-aerobic  chemistry  even  with  reduced  species  present  in  such 
low  concentrations.  Clearly,  in  a  system  experiencing  severe  anoxic 
conditions,  the  improved  simulations  resulting  from  the  inclusion  of  the 
anaerobic  subroutines  would  be  more  apparent.  Additionally,  the  ability 
to  simulate  the  presence  of  reduced  species  may  be  of  prime  importance. 
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PART  V:  SUMMARY 


78.  This  report  documents  early  testing  and  development  of  13  new 
algorithms  added  to  CE-QUAL-R1  to  simulate  anaerobic  biogeochemical  pro¬ 
cesses  in  reservoirs.  In  light  of  this  initial  evaluation,  further 
improvements  and  modifications  have  been  made  in  CE-QUAL-R1.  Some  of 
these  refinements  are  cited  here.  However,  as  the  refinements  postdate 
this  study,  they  are  not  incorporated  in  it.  Thus,  this  report  also 
serves  as  an  introduction  to  the  algorithms  and  as  a  guide  to  calibrat¬ 
ing  the  model  for  further  studies. 

79.  Two  important  simplifying  assumptions  based  on  field  data  and 
other  similar  modeling  studies  facilitate  predictions  related  to  water 
quality  under  low  DO  conditions.  First,  a  DO  concentration  of  0.5  mg/Jl 
in  the  water  column  is  used  to  initiate  these  processes.  Second, 
products  of  anaerobic  processes  are  released  simultaneously  from  sedi¬ 
ments  using  zero-order  rate  kinetics. 

80.  In  order  to  develop  good  predictions  of  the  effects  of 
anaerobic  conditions,  temperature  and  DO  predictions  must  first  be 
accurate.  Prior  to  calibration  and  testing  of  the  portions  of  the  model 
dealing  with  anaerobic  conditions,  thorough  examinations  of  the  vari¬ 
ables  affecting  thermal  structure  and  depletion  of  DO  in  the  water 
column  were  undertaken.  In  this  series  of  calibrations,  the  DO  pre¬ 
dictions  for  the  model  of  DeGray  Lake  proved  very  sensitive  to  changes 
in  DOM  decay  rate  and  somewhat  less  sensitive  to  detritus  and  sediment 
decay  rates. 

81.  Numerous  simulations  were  made  to  calibrate  CE-QUAL-R1  with 
the  new  algorithms  using  1979  data  from  DeGray  Lake  field  studies. 
Overall,  concentrations  of  anaerobic  products  were  most  sensitive  to 
changes  in  sediment  release  rates.  Problems  with  simulations  of  ammonia 
and  phosphate  concentrations  have  been  resolved  by  inclusion  of 
adsorption  algorithms  (Wlosinski  and  Collins  1985a).  Calibrations  of 
metalimnetic  oxygen  depletion  gave  realistic  results;  hypolimnetic  pro¬ 
files  deviated  more  from  field  observations.  Some  differences  can  be 
reconciled  by  recognizing  that  field  data  for  comparison  were  obtained 


from  a  shallower  site,  farther  up  the  reservoir.  Most  of  the 
hypolimnetic  deviations  can  be  interpreted  as  consequences  of  relatively 
weak,  benthic,  anaerobic  conditions,  typical  of  DeGray  Lake,  that  create 
low  concentrations  of  reduced  species  localized  near  the  bottom  which 
are  difficult  to  resolve  with  the  model.  Recent  testing  of  the  model  on 
Eau  Galle  Lake,  Wise.,  where  anaerobic  conditions  were  more  severe, 
produced  better  comparisons  with  less  calibration  (Wlosinski  and 
Collins  1985b) .  In  spite  of  low  concentrations  of  reduced  materials  in 
DeGray  Lake,  the  model  reflects  their  slight  impact  on  DO 
concentrations. 

82.  Discrepancies  between  modeled  and  observed  conditions  in  the 
confirmation  simulations  for  1980  can  be  ascribed  to  the  same  causes  as 
in  the  calibration  work.  Differences  in  conditions  early  in  1980  also 
contributed  to  predicted  water  quality  deviating  from  actual  system 
behavior.  In  particular,  in  early  February  1980  the  reservoir  was 
isothermal  but  still  chemically  stratified.  However,  the  isothermal 
condition  caused  the  model  to  Induce  mixing  prior  to  its  occurrence  as 
indicated  by  field  data.  This  result  indicated  the  need  to  examine 
criteria  controlling  mixing.  Adjustments  of  the  relevant  coefficients 
might  have  improved  results  for  1980  without  compromising  1979  results. 
Overall,  trends  and  concentration  peaks  match  field  data  well,  particu¬ 
larly  in  the  metalimnion  which  is  more  important  for  the  DeGray  Lake 
system;  these  results  can  be  attributed  to  calibrating  the  model  for 
metalimnetic  conditions  which  become  anaerobic  earlier  and  more 
severely. 
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